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Abstract 

The combined effect of mean flow and rotation on hexagonal patterns is investigated 
using Ginzburg-Landau equations that include nonlinear gradient terms as well 
as the nonlocal coupling provided by the mean flow. Long- wave and short-wave 
side-band instabilities are determined. Due to the nonlinear gradient terms and 
enhanced by the mean flow, the penta-hepta defects can become unstable to the 
induced nucleation of dislocations in the defect-free amplitude, which can lead to 
the proliferation of penta-hepta defects and persistent spatio-temporal chaos. For 
individual penta-hepta defects the nonlinear gradient terms enhance climbing or 
gliding motion, depending on whether they break the chiral symmetry or not. 
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1 Introduction 



Patterns in Rayleigh-Benard convection have been extensively explored, both 
theoretically and experimentally, as paradigms to investigate spatio-temporal 
chaos and transitions from ordered to disordered states. Arguably, one of the 
most interesting chaotic states is that of spiral defect chaos. For small Prandtl 
numbers, it arises from roll convection at moderate heating rates in systems 
with large aspect ratio [1] . In this complex chaotic state spirals, disclinations, 
and dislocations are persistently created and annihilated. It is found to be 
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driven by a large-scale mean flow, which arises due to the curvature of the 
convection rolls and becomes prominent for small Prandtl numbers. 

Another classic chaotic state that arises from a planform of convection rolls 
is the domain chaos resulting from the Kiippers-Lortz instability in rotating 
convection [2-4]. In this instability rolls are unstable to rolls with a similar 
wavelength but different orientation. The resulting state is characterized by 
domains of convection rolls with different orientation that persistently invade 
each other. 

Motivated by the strong impact that mean flow and rotation have on the 
stability of rolls and by the interesting chaotic states that result from it, we 
consider in this paper the effect of mean flow and rotation on the stability 
and dynamics of hexagon patterns and their defects. Hexagon patterns are 
common phenomena in various pattern-forming systems and are readily ob- 
tained also in Benard-Marangoni convection driven by surface tension (e.g. 
[5,6]) and in buoyancy-driven non-Boussinesq convection (e.g. [7]). In both of 
these convection systems mean flow is expected to be important for small and 
moderate Prandtl numbers. 

In the absence of mean flows, the side-band instabilities and dynamics of 
hexagons with broken chiral symmetry (e.g. by rotation) have been investi- 
gated in the context of three coupled Ginzburg-Landau equations [8] as well 
as in a long-wave model for Marangoni convection [9] and a model of Swift- 
Hohenberg type [10]. Even on the level of the lowest-order Ginzburg-Landau 
description rotation makes the system non-variational, and oscillatory ampli- 
tude [11-13] and side-band instabilities appear [8]. If the nonlinear gradient 
terms are retained in the Ginzburg-Landau equations the latter can lead to 
an interesting state of spatio-temporal chaos [8]. 

In the absence of rotation, mean flow couples differently to the two steady, 
long- wave phase modes of the weakly nonlinear hexagon patterns [14]; only 
the stability limits due to the transverse phase instability are affected by the 
mean flow, while those corresponding to the longitudinal phase mode are un- 
changed. As a result, for sufficiently small Prandtl numbers the stability limit 
for large wavenumbers is given by the transverse phase mode, while for small 
wavenumbers the longitudinal mode becomes dominant. This is particularly 
interesting, since without mean flow the longitudinal mode is usually of little 
importance [15,14]. While the transient patterns arising from the longitudinal 
instability remain quite ordered, those ensuing from the instability involving 
the transverse mode typically exhibit domains of hexagons with quite different 
orientation [14]. 

In this paper, we investigate the combined effects of rotation and of mean 
flow on weakly nonlinear hexagonal patterns using appropriately extended 
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Ginzburg-Landau equations. We retain all three possible cubic non-linear gra- 
dient terms. The nonlinear evolution of the side-band instabilities leads to 
the formation of dislocations that later combine to penta-hepta defects. In 
parameter regimes in which the nonlinear gradient terms, rotation, and mean 
flow are significant we find quite intriguing defect dynamics: the presence of 
penta-hepta defects induces the nucleation of a dislocation pair in the defect- 
free amplitude, which can lead to the proliferation of defects and to persistent 
spatio-temporal chaos. We further study the effect of the nonlinear gradient 
terms on the motion of single penta-hepta defects by calculating their mobility 
and the Peach-Kohler-type force acting on them. 

This paper is structured as follows. In Section 2 we formulate the problem 
by extending previous work on mean flow in hexagons [14] to include the 
breaking of the chiral symmetry by rotation. Then we investigate in Section 
3 the linear stability of hexagons with respect to side-band perturbations. We 
demonstrate the induced defect nucleation and the resulting spatio-temporal 
chaos in Section 4. In Section 5 we investigate the effects of the nonlinear 
gradient terms on the motion of isolated penta-hepta defects both analytically 
and numerically. We discuss our findings in Section 6. 



2 Amplitude equations 

For small Prandtl numbers and with rotation the usual Ginzburg-Landau 
equations for the three complex amplitudes Aj making up the hexagon pat- 
terns need to be extended in two ways. The mean flow, which is driven by 
long-wave modulations of the convective amplitude, provides a nonlocal cou- 
pling of the roll modes [16]. Rotation breaks the chiral symmetry. This is 
reflected in the difference between the cubic term coupling A\ to A2 and that 
coupling A\ to A3 [11,12,8]. Furthermore, to cubic order it introduces an ad- 
ditional nonlinear gradient term in the equation for the amplitudes and an 
additional term in the leading-order equation for the mean flow Q [17]. The 
equations for the three hexagon modes in rotating convection at finite Prandtl 
numbers thus read 



d t A 3 = (/i + (hj • V) 2 - \Aj\ 2 -{v- 7)|^Vi| 2 - (u + -y)\A j+1 1 2 ) Aj 
+A*_ 1 A* j+1 -if3A j (T j -V)Q 

- a 3 )A*_ 1 (h j+1 ■ V)A* +1 + i( ai + a 3 )A* j+1 (n j _ 1 ■ V)A*_ 



+ia 2 (A*_ 1 (r j+1 ■ V)A* +1 - A* +1 {r^ ■ V)A*_ X ), 
V 2 Q = ]T [2(ni ■ V)(fi • V) + r (ifn ■ V) 2 - (n ■ V) 2 )] \A t 




(2) 



i=i 
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with j = 1, 2, 3 cyclically permuted in (1) and fij and fj denoting unit vectors 
parallel and perpendicular to the critical wave-vector associated with ampli- 
tude Aj, respectively. With respect to the rotation rate, the coefficients 7, 
o?3, and r are odd functions, while the other coefficients are even. They also 
depend on other physical parameters such as the Prandtl number. Note that 
in the presence of rotation not only the difference ((n^ • V) 2 — (fj ■ V) 2 )|Aj| 2 
but also the corresponding sum is allowed by symmetry, which corresponds to 
V 2 |Aj| 2 and leads to a contribution Qi oc X)f=i |A| 2 to Q. Upon insertion in 
(1) Qi contributes to the cubic nonlinear gradient terms and provides there- 
fore a local rather than a non-local coupling of the roll modes. In this paper 
we focus on the non-local coupling and neglect this term along with the other 
cubic gradient terms. 



The above equations allow for several stationary, spatially periodic patterns. 
Rolls with amplitudes Aj = y/Jl — q^ e iqnj ' x , A j±1 = 0, exist for ji > q 2 , and 
are stable to homogeneous perturbations for \i > fiR + q 2 , where /ir — (1 + 
2qa\) 2 /(l — v) 2 . Hexagon solutions with wavenumber slightly different than 
the critical value q c are given by Aj = R e tqn: >'* with 



(1 + 2q ai ) ± ^(1 + 2qa 1 ) 2 + 4(/x - q 2 )(l + 2v) 

R ° ~ 2(TT2^ • (3) 



We will consider only such equilateral hexagons, for which the wavenumbers 
in all three modes are equal. Non-equilateral hexagon patterns have been 
discussed in the absence of mean flow and rotation in [18]. Since the quadratic 
coupling coefficient in equation (1) has been scaled to +1, the hexagon solution 
corresponding to the amplitude R with a minus sign in front of the square 
root is always unstable. In the following we will consider only the hexagon 
solution of amplitude Rq with the plus sign in equation (3). It is stable to 
homogeneous perturbations if both u = Rq(1 — u) + (1 + 2aiq)R > and v = 
2R 2 ) (l + 2v)-(l + 2a 1 q)R > 0. Mixed modes with amplitudes Aj = l/(v-l) 
and Aj±i = ^J(fi — q 2 — A 2 ) /(l + v) also exist but are always unstable with 
respect to rolls or hexagons. The stability of these stationary, spatially periodic 
solutions with respect to homogeneous amplitude perturbations is summarized 
in the bifurcation diagram in Fig.l: the hexagons first appear at v — (or 
\x — Hsn in the bifurcation diagram) through a saddle- node bifurcation, and 
become unstable to 'oscillating hexagons' via a Hopf bifurcation at u — (or 
fi = in Fig.l) with a Hopf frequency to c = 2v^37(l + 2qa\) 2 / {y — l) 2 . Here 
we focus on the instability of steady hexagons in the range fi sn < /j, < fin', 
a detailed discussion of the oscillating hexagons for ji > fin can be found in 
[11,12,8,19]. 
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Fig. 1. Sketch of bifurcation diagram of simple, spatially periodic solutions of (1,2) 
with rotation. 

3 Phase equation and general stability analysis 



Following the procedures in [20,8], we first derive the phase equation that 
describes the long-wave side-band instabilities of hexagons. The slightly per- 
turbed hexagon solution 

A, = Roe^ x (l + r 3 + i^), j = 1, 2, 3 (4) 



is substituted into equations (1,2), where rj and <pj are the small amplitude 
and phase perturbations, respectively. We introduce the super-slow scales d t — > 
8 2 d t and V — > 5V with \8\ <C 1, and adiabatically eliminate the perturbations 
in the amplitudes and in the total phase $ = 0i + 2 + 3 by expressing 
them in terms of the two translation phase modes <f) x = — (02 + 03) an d 
<p y = (02 — 03)/y / 3. At linear order in 5, the mean-flow amplitude Q can be 
written in terms of the phase vector = (<f> x , <p y ) as 



Q = d x2 V ■ + d ± e z ■ V x 0, (5) 
where 



d± = 4^32 K« - ^ w + (™ + |)>/3asi2g], (6) 

dx2 = 4 U 2 + ^2 + 2 ) w ~ ( M ~ 2 r)V3a 3 Rl], (7) 

iu = 2gi2g + (an + A/3a 2 )i2o, (8) 

^ = 2^07, (9) 

and e z is the unit vector perpendicular to the x — y plane (following a right- 

— # 

hand rule). The linearized phase equation for thus reads 



dt$= D° V 2 + (D? { - D° ± )V(V • 0) - L>° xl (e, x V 2 0) + £>£ 2 (e a x V)(V0) 
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-/3V x (Qe z ), 



(10) 



where the coefficients with superscript correspond to the infinite Prandtl 
number case (/? = 0) [8]: 



D u ,=- + 



1 1 \Rlu 



o< 



4 4u 2 + uj 2 2 



(ai + v^3«2) 2 + 3a|] — V3RoLua>3q — 2-ug 2 j , 



=D° ± + I - -{Rla^cn - >/3a 2 ) - -Ro(3«i - + 2g 2 } , 

-°xi = -I - 5 9 \ T^-^o ( a i + v / 3«2) 2 + 3ao + 2v / 3-R wa3<? - uq 2 \ , 

^°x2 = (>/3i^ai - VSRoq) ■ (11) 
Substituting equation (5) into the linear phase equation (10), we obtain 

d t $= D ± V 2 + (£>u - D ± )V(V ■ $) - D xl (e z x V 2 0) + ^x 2 (e, x V)(V0),(12) 

where D ± = D° ± + (3d ± , D\\ = £>}, D x i = £>°i and £>x 2 = £>°2 + /^ x2 - 
The terms (3d± and /3d x2 are the mean flow contributions to the diffusion 
coefficients. 



The two eigenvalues for normal mode solutions to equation (12) are 

1.2 r 



Dm + D ± ± - D ± ) 2 - 4D xl (D xl - D 



x2) 



(13) 



where k is the magnitude of the wave- vector of the normal-mode perturbations 
proportional to e at+lk ' yi . The phase instability becomes oscillatory when the 
discriminant 

_ = (£, _ D±) 2 _ 4D xl ( Dxl _ Dx2 ) ( M ) 



is negative, which is possible only when the chiral symmetry is broken. 

The stability of the hexagonal pattern to general (including short-wave) per- 
turbations is examined by solving linearized equations similar to those in 
[14]. Without rotation and the nonlinear gradient terms, the general stabil- 
ity boundaries correspond to a steady bifurcation (real eigenvalues), and they 
always coincide with the long- wave stability boundaries [13,14]. This is no 
longer true if rotation and nonlinear gradient terms are included. Hexagons 
may undergo instability via short-wave instabilities, which may be steady or 
oscillatory depending on the parameters [13]. In the following stability dia- 
grams, we display both the stability boundary for the long-wave perturbations 
(lines) and the short-wave stability boundaries (circles). 
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Wavenumber q Wavenumber q Wavenumber q 



Fig. 2. Stability diagrams for v = 2, a\ = a 2 = «3 = 0, 7 = 0.05. a) (3 = r = 0, 
b) (3 = — 1 and r = 0, c) (3 = — 1 and r = 1. Circles denote short-wave instabilities. 
Thick lines denote long- wave instabilities: oscillatory (solid) and steady (dashed). 
Dotted line denotes saddle-node bifurcation, thin solid line gives Hopf bifurcation 
to oscillating hexagons. 

Figs.2a,b,c represent stability diagrams for different values of (3 and r with 
7 = 0.05 and «i = a 2 = a 3 = 0. Fig.2a is for (3 = and r = 0, Fig.2b for (3 = 
— 1 and t = 0, and Fig. 2c for /3 = — 1 and r = 1. In all the following stability 
diagrams the thick dashed curves (labeled 1, 3 and 4 in Fig. 2) correspond to 
steady phase instabilities, whereas the thick solid lines (labeled 2) denote the 
oscillatory phase instability. The symbols denote short-wave instabilities. As 
is the case without rotation [14], mean flow renders the stability boundaries 
asymmetric with respect to the band-center (q = 0) . In addition, we find that 
the 'bubble' enclosed by curve 1 expands as the mean flow becomes larger in 
amplitude (figs. 2a and 2b), suggesting that the mean flow tends to diminish 
the importance of the oscillatory modes. For (3 < 0, this effect appears to 
be more prominent for negative q than for positive q and can eliminate the 
oscillatory instability altogether (cf. Figs. 3b and 4b below). 




Wavenumber q Wavenumber q 

Fig. 3. Stability boundaries for 7 = 0, = 0, 03 = 0.7, v = 2. a) /3 = 0, b) 
(3 = —3. Circles denote steady short-wave instability. Solid thick line gives oscillatory 
long-wave instability, dashed thick line steady long-wave instability. 

In the stability diagrams shown in Figs.3a,b, we focus on the nonlinear gradient 
term that is due to rotation (03 = 0.7) and set the other nonlinear gradient 
terms as well as the other rotation terms to zero (7 = 0, r = 0). While the a\- 
term makes the stability limit of the hexagons to the mixed-mode asymmetric 
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Wavenumber q Wavenumber q 

Fig. 4. Stability boundaries for v = 2, a± = «2 = 0, «3 = 0.7, 7 = 0.2, r = 0.5. a) 
(3 = —0.1, b) f3 = —3. Circles denote steady short-wave instability. Solid thick line 
gives oscillatory long-wave instability, dashed thick line steady long-wave instability. 



in q and can shift the transition from hexagons to rolls to large values of /i [21], 
the term corresponding to 0:3 does not affect that amplitude instability. For 
(3 = (Fig. 3a) the long-wave instability is oscillatory, while the short-wave 
instability is steady. When [3 is decreased to [3 = — 3 (Fig. 3b) the oscillatory 
instability is completely suppressed for q < and replaced by the two steady 
phase instabilities (labeled 3 and 4), whereas for q > it is still there (2), but 
it is mostly preempted by a steady short-wave instability. Only for very small 
values of /1 the oscillatory instability remains relevant. 

If in addition to a.3 also the cubic rotation term 7 is present the stability 
diagram becomes asymmetric even in the absence of mean flow. This is illus- 
trated in Fig. 4a, which gives the stability limits for 7 = 0.2, a% = 0.7 and 
t = 0.5, with a\ = = ct2 and (3 = —0.1. Note that for small (3 the stability 
boundaries depend very little on (3 and are indistinguishable on the scale of 
Fig. 4a for (3 in the range —0.2 < (3 < 0. As will be discussed in section 4, 
however, the nonlinear evolution of the pattern due to the instabilities can 
depend substantially on f3 even in this regime. The origin of the asymmetry 
in q can be seen from the diffusion coefficients and the growth rate given in 
(11) and (13), respectively. For 07 = = «2 the only terms that are odd in q 
involve the product a^uo. For these parameter values the short-wave instability 
is steady and the long-wave instability is oscillatory. The steady long-wave in- 
stability (short segment of a dashed line) has been shifted up in /i all the way 
to the amplitude instability to oscillating hexagons and is preempted by the 
long-wave oscillatory instability. In Fig.4b the mean flow strength is increased 
to (3 = —3. Then the stability limits are given solely by steady long-wave 
instabilities and hexagons are stable only over a small range in ji. 



8 




Fig. 5. Temporal evolution of the number of dislocations for different values of [3. 
Other parameters as in Fig. 4a. (a) j3 = for system size 42 x 48 (dashed line) and 
84 x 96 (solid line shows quarter of the number of dislocations) (b) ft = —0.1 (dashed 
line) and (3 = —0.2 (solid line) for system size 42 x 48. 

4 Defect Proliferation 



We numerically simulate equations (1,2) to investigate the non-linear dynamics 
ensuing from the linear instabilities and focus on the combined effect of the 
nonlinear gradient terms and the mean flow. Without rotation, the mean flow 
couples only to the transverse phase mode and makes it possible to discern 
its instability from that of the longitudinal phase mode [14]. Both instabilities 
lead to the formation of PHD and eventually return the wavenumber to the 
stable band. While the transient disorder generated by this process is quite 
different for the two instabilities, they both, in the absence of the nonlinear 
gradient terms, lead to an essentially monotonic decay of the defect number 
after it has reached its initial maximum. 

With the nonlinear gradient terms included the defect dynamics can become 
much more complex. In Fig.5a,b we show the temporal evolution of the number 
of dislocations for parameters corresponding to the stability diagram shown 
in Fig. 4a with /x = 1, which are characterized by a relatively strong nonlinear 
gradient term that breaks the chiral symmetry («3 = 0.7)0. Three values for 
the mean flow strength are used, (3 = 0, (3 = —0.1, and (3 = —0.2. The initial 
condition is a periodic hexagon pattern with wavenumber q = —0.6 (diamond 
in Fig. 4a) that is perturbed with a small long-wave modulation. 

Since the initial wavenumber is only slightly outside the linear stability bound- 
ary it takes quite a long time until the first defects are generated. Then the 

2 While for the parameters corresponding to fig. 4b hexagons of all wavelengths are 
unstable for fi > 0.6 numerical simulations in this regime show that the defects 
arising from the instability serve as nucleation sites for rolls, which then invade the 
whole system. 
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Fig. 6. Root-mean-square of the mean flow velocity as a function of time for 
(3 = —0.1 and j3 = —0.2. Other parameters as in Fig.4a. 

number of dislocations reaches a maximum very rapidly. Without mean flow 
(P = 0, Fig. 5a), it decays subsequently to very small values and the final 
state is presumably stationary (cf. Fig. 11 below). In contrast to the case with- 
out nonlinear gradient terms and without rotation [14], further defects are, 
however, created throughout the simulation time, in particular in the larger 
system of size 84 x 96 (solid line). 

When turning on the mean flow in addition to the nonlinear gradient terms, 
the evolution becomes more complex as shown in Fig. 5b. Decreasing (3 from 
P — to P — —0.2 shifts the stability limits very slightly towards less negative 
wavenumbers and renders the hexagons slightly more stable. The initial rise 
in the number of dislocations is therefore delayed to later times and the initial 
peak is smaller. Even though this change in the stability limits cannot be 
discerned on the scale of Fig. 4a, the weak induced mean flow is sufficient to 
lead to strong and persistent fluctuations in the number of defects (solid line 
in Fig. 5b). In fact, somewhat smaller and slower fluctuations persist even for 
ft = —0.1 (dashed line). The strong correlation between the number of defects 
and the mean-flow strength is apparent when comparing Fig. 5b with Fig.6, 
which shows the spatial average of the root-mean-square of the mean-flow 



4-1 Induced Defect Nucleation 

A closer inspection of the defect dynamics shown in Fig. 5 reveals that the 
temporal fluctuations result from the induced nucleation of new dislocations 
in the vicinity of pre-existing penta-hepta defects. One such event is shown 
in Fig. 7. The initial configuration consists of a PHD denoted by (0, — ,+), 
which is comprised of a dislocation with negative topological charge in mode 
A 2 and one with positive charge in mode A 3 (grey circles in Fig. 7). Mode Ai 



velocity J< (VQ) 2 > as a function of time. 
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(a)Time = 881.25 (b)Time = 887.5 




(c)Time = 893.75 (d)Time = 900 




Fig. 7. Induced nucleation of dislocations by penta-hepta defects for (3 = —0.2 and 
[1 = 1. Other parameters as in Fig. 4a. Solid lines are zero contour lines for real part 
of A\, and dashed lines for the imaginary part. Grey circles are dislocations in A<i 
and A% and the arrows indicate the mean flow velocity field. 

does not have any dislocations in the vicinity of this PHD. It is, however, 
perturbed due to the presence of the PHD and the zero-contour lines of the 
real and the imaginary part of Ai (thick solid and dashed lines) are twisted 
around the dislocations in A2 and A3 (Fig. 7a). This twisting reflects to some 
extent the advection of this mode by the mean flow (thin lines). Soon two 
dislocations appear in A\ (black circles in Fig. 7c), which then bind with the 
dislocations constituting the initial PHD to form two new PHDs, (+.— , 0) and 
(— , 0, +). The new PHDs typically move apart and can induce the nucleation 
of further dislocations leading to the proliferation or multiplication of defects 
as sketched in Fig.8. There the further nucleation of a pair of dislocations in 
A 2 and A 3 leads eventually to four PHDs. Such proliferation processes have 
recently been found also in simulations of Kuramoto-Sivashinsky equation and 
of Ginzburg-Landau equations without rotation and without mean flow but 
with the nonlinear gradient terms corresponding to cti.2 [22]. 

The nonlinear gradient terms appear to be central for the induced nucleation 
in that they lead to a separation of the dislocations making up the penta-hepta 
defects as is illustrated in Fig. 7a. The dependence of the distance between the 
dislocations on 03 is shown in Fig. 9. For a 3 > 0.6 the PHD becomes unstable 
through the induced nucleation of dislocations. We find that the dislocations 
are also separated if only the nonlinear gradient terms corresponding to oti 
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Fig. 8. Sketch of defect proliferation by induced nucleation. 
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Fig. 9. Distance between dislocations within a bound PHD as a function of 03 for 
(3 = 0. Other parameters as in Fig. 4. 

or a 2 are included (see also [22]). Breaking the chiral symmetry through the 
7-term does not affect the distance between the dislocations in a PHD pair, 
nor does the mean flow. 




The stability of the PHD, however, depends not only on the distance between 
its two constituent dislocations but also on the mean flow. This is shown in 
Fig. 10 where the stability limit of PHDs is given as a function of a 3 and j3 for a 
background wavenumber q — 0. It is obtained by direct numerical simulations 
of a PHD pair in a system of size 42 x 48. Note that for this system size 
there is still a small, but noticeable interaction between the two PHDs that 
have to be placed in the system to satisfy the periodic boundary conditions. 
In the stable regime the two PHDs move relatively slowly towards each other 
and would eventually annihilate each other, whereas in the unstable regime 
the seeded PHDs nucleate new dislocations and get transformed into different 
PHDs before they reach each other. 

Clearly, increased mean flow greatly facilitates the induced nucleation of dis- 
locations rendering the PHDs much less stable even though the distance of the 
constituent dislocations is hardly affected by the mean flow. Thus, even rela- 
tively small values of a% can be sufficient to induce nucleation if \/3\ = 0(1). 
When «3 is decreased below 013 0.1 the induced nucleation as described 
above is replaced by a different process in which new dislocations appear not 
in the previously defect-free mode but in the modes that carry already a dis- 
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a 3 



Fig. 10. Stability limits of penta-hepta defects (squares) and limit of persistent 
spatio-temporal chaos (diamonds) as a function of mean-flow strength and of 013. 
Remaining parameters are fx = 1, v = 2, 7 = 0.2, r = 0.5, and q = 0. 

location. This is presumably due to the fact that with increasing mean-flow 
strength the side-band stability limit of the periodic hexagon pattern is shifted 
further to the left (for example, see Fig.l in [14]) and comes very close to the 
background wavenumber q = employed in Fig. 10. 



4.2 Persistent PHD- Chaos 



A particularly interesting aspect of the simulations shown in Fig. 5b is the fact 
that the induced nucleation of dislocations is not merely transient as had been 
the case in [22], but instead persists for the whole duration of the simulation. 
To bring out the persistence of the dynamics more clearly, we give a detailed 
analysis of the pattern evolution by measuring the local wavenumber q^ (x) of 
each component, 



^eRI^I for j = 1,2, 3. (15) 



We extract the spatial average of the transverse and of the longitudinal com- 
ponents qij and q t j of q^, respectively, 



Sr % ' d x _ Jr f? ' d x 



<to = Vq^ /rd2x , <* = *r*= /rd2x • (16) 

Here V denotes the spatial domain. Figs.lla,b show the temporal evolution of 
the two wave-vector components corresponding to the evolution of the num- 
ber of defects shown in Figs.5a,b (for system size 42 x 48). The rapid ini- 
tial change in both components corresponds to the period in time when the 
long-wave instability first generates dislocations. For all three values of (3 the 
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Fig. 11. Spatially averaged local wave vectors as functions of time: thick lines are 
the transverse components qtj and thin lines are the longitudinal components qjj of 
the wave- numbers. Solid lines are for A\, dashed lines for Ai and dotted lines for 
A3. Panel (a) is for (3 = and panel (b) is for (3 = —0.1 and (3 = —0.2. 

longitudinal component rapidly relaxes to a very small value near qu = and, 
with the exception of small fluctuations, stays there for the whole duration 
of the simulations. Similarly, for (3 = the transverse component reaches in 
a somewhat slower transient a stationary value of q[J ~ 1.6. The resulting 
reconstructed pattern at the final time t = 2500 is shown in Fig. 12a. It still 
contains a few PHD, which move very slowly. Compared to the initial pattern, 
in which the hexagons were aligned with the x-axis, the pattern is rotated by 
a small amount reflecting the change in the average wave- vector, in particular 
its non-zero transverse component. 

For (5 = —0.1 and (5 = —0.2 the number of PHDs keeps fluctuating throughout 
the simulation. The corresponding evolution of the average wave-vector is 
shown in Fig. lib. As is the case for (3 = 0, during the initial phase of defect 
generation the longitudinal component relaxes rapidly to qfj = while the 
transverse component reaches a value of qtj ~ 1 during that phase. Thereafter, 
the longitudinal component remains near whereas the transverse component 
keeps growing at a rate which increases with decreasing (negative) (3. In the 
reconstructed pattern, shown in Fig. 12b at t — 2500, this manifests itself 
in a reduced overall wavelength of the hexagons and, more importantly, an 
increased rotation of the pattern. Following the reconstructed pattern over 
time one can see that on average it is steadily rotating in a counter-clockwise 
fashion. 

To interpret the evolution shown in Fig. lib it is useful to write Aj(x, t) as 



For constant Aj this corresponds to a hexagon with non-vanishing transverse 
wave- vector components that are equal in all three modes. Insertion of (17) 



Aj = Aj(x,t)e 



,iqtfj-x 



J = 1,2,3. 



(17) 
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Fig. 12. Reconstructed hexagonal patterns at the end of the simulations (t = 2500) 
for (3 = on the left and (3 = —0.1 on the right. 

into (1,2) shows that the solution and its stability does not depend on q t . In 
particular, the contributions from «2 cancel each other. This independence 
of q t reflects the isotropy of the underlying physical system and the leading- 
order approximation of the critical circle by its tangents in the direction of fj 
at each of the three critical wave- vectors corresponding to the modes A y This 
suggests that the steady increase of eft] with fixed qjj should be interpreted as 
the representation of a rotation of the physical pattern (at fixed magnitude of 
the wave- vector) within the approximation of the Ginzburg-Landau equations 
(1,2). A similar phenomenon had been observed for hexagons with broken 
chiral symmetry in the absence of mean flow [8]. In that system it had been 
found that if the operator fij ■ V was replaced by the Newell- Whitehead-Segel 
operator [23,24] or by the Gunaratne operator [25] the steady growth of the 
wave-vector components did not follow fj but instead followed the respective 
lines in Fourier space along which the growth rate of perturbations of the basic 
state is constant within these two different approximations [26]. 

Given the above discussion, it is appropriate to consider within the approx- 
imation (1,2) all hexagon solutions with gy = as having a wave-vector at 
the band-center, with q t j indicating the orientation of the pattern in space. 
This raises the question why new dislocations are generated persistently even 
though the background wavenumber is in the band center where the periodic 
hexagon patterns are linearly stable? In [22] it had been found that without 
mean flow the PHDs can be unstable even if the wavenumber of the back- 
ground hexagon pattern is inside the stable band, but away from the band 
center. As was shown in Fig. 10 above, in the presence of mean flow PHDs 
can be unstable even at the band center. We have not investigated whether 
there are parameter regimes for which the PHDs are unstable at all back- 
ground wavenumbers. While instability at the band center is consistent with 
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the persistence of spatio-temporal chaos since the background wavenumber of 
the chaotic state is q = 0, it is not sufficient. This is indicated in Fig. 10, where 
also the persistence limit for the chaotic state is shown. Below the squares the 
PHD's are unstable but chaos only persists for values of (3 below the diamonds. 
In the parameter regime between these two lines defects are being created in- 
termittently during the transient towards the stationary state, which leads to 
fluctuations in the defect number, but the defect creation rate is too small 
compared to the annihilation rate to sustain the chaos. 

To illustrate explicitly how the instability of a PHD at the band center can lead 
to persistent PHD chaos, we performed simulations that start from hexagons 
at the band center (q = 0) with a pair of PHDs as a seed for the nucle- 
ation process. The resulting evolution of the number of defects and of the 
spatially averaged transverse wave- vector component is shown in Figs.l3a,b, 
respectively. Compared to the simulations shown in Fig. 12b, the value of a 3 is 
reduced to 0:3 = 0.5. As a consequence even for f3 = —0.2 no new dislocations 
are generated. For (3 = —0.4 there is an initial volley of induced nucleation, 
but eventually the defects annihilate each other again completely. Only for 
(3 < —0.5 indications for persistent nucleation are seen. This illustrates that 
the instability of PHD is not a sufficient condition for persistent chaos; it is 
necessary that creation balances annihilation for some finite average number 
of PHDs. 

Worth noting is the precipitous decline in the defect number for (3 = —0.6 
near t = 900. At that time the transverse wave-vector component reaches 
values close to the maximal spatial resolution for the number of modes used 
(128 x 128). More detailed studies of the case with a = 0.7 and (3 = —0.2 
have shown that if the number of modes is increased the suppression of the 
induced nucleation is delayed until larger values of gTJ are reached. Thus, the 
saturation of ~qTj for (3 = —0.5 in Fig. 13b and the associated end of the chaotic 
activity is a numerical artifact. Note, however, that the chaotic defect states 
reported here have lasted substantially longer than the transient states found 
in [22]. 

Since the defect nucleation persists even for wavenumbers near the band-center 
it is not surprising that similar chaotic states are reached when the stability 
limit for periodic hexagons is crossed on the large-g side. It should be noted 
that in simulations that started from a slightly perturbed hexagon pattern 
beyond the low-g stability limit no persistent nucleation was found when the 
wavenumber was too far in the unstable regime. In that case only a very 
large number of defects was generated in the initial phase, which then quickly 
annihilated each other. 

Of course, the steady increase of the transverse wave- vector component appar- 
ent in Fig. lib implies that after a finite amount of time the magnitude of the 
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Fig. 13. a) Temporal evolution of number of dislocations for a = 0.5 and various 
values of (3. Other parameters as in Fig.5. b) Evolution of the spatially averaged 
transverse wave vector component qij. 

transverse component qtj becomes of the same order as the critical wavenum- 
ber. At that point the Ginzburg-Landau equations clearly will have become 
invalid, since they are based on the approximation that the amplitudes A,- 
vary only slowly compared to the critical wavelength. We have therefore also 
simulated a modified Swift-Hohenberg equation that corresponds to the same 
amplitude equations with coefficients corresponding to the values used in the 
simulations presented here. As expected, we find that the hexagon pattern 
keeps rotating on average as the defect proliferation and annihilation continue 
indefinitely [27]. 

In general, the other nonlinear gradient terms corresponding to a,\ and «2 have 
to be taken into account as well. We have performed simulations for selected 
parameter values of a\ and a 2 that correspond, for instance, to realistic values 
for surface-tension-driven Benard-Marangoni convection («i = 0.3 and a.2 = 
1.1, corresponding to a Prandtl number ~ 10 in a single-fluid model [28]; 
«i = —0.1 and «2 = 0.4 corresponding to a two-fluid model a water-air layer 
[29]). While we find induced nucleation even in the absence of rotation and 
without a 3 , persistent spatio-temporal chaos arises only if a 3 > 0.8 in both 
cases (f3 = -0.2). 



5 Effects of non-linear gradient terms on PHD motion 

The motion of individual PHD's within the lowest-order Ginzburg-Landau 
equations has been studied in detail in [30,31], where in extension of the results 
for dislocations in roll patterns [32] the dependence of the velocity of a PHD on 
the background wave- vectors of the three hexagon modes has been determined 
semi-analytically. The effect of the mean flow on the defect motion has been 
discussed in some detail in [14]. Here we turn to the impact of the other non- 
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variational terms, specifically the non-linear gradient terms, on the motion of 
the PHD's. 



In our analysis we closely follow the approach of [31]. We factor out the back- 
ground wave- vectors of the three hexagon modes by writing Aj = a^e"^'*. 
Assuming a constant defect velocity V, the time derivative in (1) is replaced 
by —V • V. We then project (1) onto the two translation modes, i.e. for each 
j we multiply (1) by d x a* and d y a*, respectively, and add all three equations 
and their complex conjugates. The resulting projection can be written as 

^ra^c IxyVy Pi ^yx^x lyy^y G, (IS) 

where the components of the mobility tensor are given by 




and < • • ■ >= // • • • dxdy denotes the integral over the domain. The terms 
on the right- hand- side of equation (1) contribute to F and G with 



F = _2 7 (S H^k+il 2 ) -i(^2qid x a*(ni ■ V)a, - c.c. 



\i=i 



—i- 



d x a*ai(fi ■ V)Q - c.c. 



\i=l 
/ 3 



2 d X 



(^2 d * a i [ a *-i(n+i ■ V)a* +1 - a* +1 (fi_i • V)a-_! 



\i=i 

I 3 



■ V)a; 



+i + <+i(irV)<-i 



— c.c. 



— c.c. 



a*_i(n i+1 ■ V)a* +1 - a* +1 (r^_i ■ Vja*^ 



c.c. 



= 7 F 7 + q F q + pFfs + aiFi + a 2 F 2 + a 3 F 3 . 



(20) 



G is defined analogously to F with d x a* in front of the square brackets in each 
of the integrals replaced by d y a*. The projection qi of the wave- vector q onto 
hi is given by q t — q^ ■ hi. To leading order the transverse components • of 
the wave- vectors do not affect the velocity of the defect. Here we only focus on 
situations where qi = go is the same for all three amplitudes. Note, that the 
other terms in (1) do not contribute to (20), which can be seen by integrating 
by parts and using the appropriate boundary conditions. 
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In the absence of the non-potential terms (7 = (3 = ot\ = a 2 = 0:3 = 0) 
the PHD is stationary for q = and the two dislocations making up the 
PHD are located at the same position. Without loss of generality, we assume 
that the PHD consists of a pair of dislocations in A 2 and A 3 . Treating the 
non-potential terms perturbatively it is sufficient to use the potential solution 
for the stationary defect to evaluate the inhomogeneous terms F and G. Far 
away from the defect core the pattern is well described by the phase equations; 
rewriting the demodulated amplitude as aj = pj(r, <f))e x9i ^\ where r is the 
radial distance from the defect core and is the azimuthal angle around the 
core, one obtains then [33] p«(r — > 00, 0) — > p and 



*i = 2^(1 -«>s(20), 



60 



1 



2V3 



e 3 = -(f> 



2V3 



1 . 2tt n 

_ + cos ( 2 - — ) 

_ + cos ( 2 + — ) 



(21) 
(22) 
(23) 



This solution is utilized for the far-field contribution to the integrals for F 
and G. The numerical solution of the Ginzburg-Landau equation (1,2) is used 
for the evaluation of the integral in the core region. It can be shown that 
the far-field contribution to the integrals G±, G2, and F 3 is zero, while that 
to the integrals iq, F 2 , and G 3 is non-zero. For example, for G\ the angular 
integral involving 9\ vanishes, and those involving 9 2 and #3 cancel each other, 
leading to a vanishing contribution to G\ from the far field. Numerically, we 
also found that the core integrals for G±, G 2 , and F 3 vanish; they are at least 
100 times smaller than the core integrals for iq, F 2 , and G 3 . In the following 
we ignore F 1 and G 7 ; similar to the integrals G±, G 2 , and F 3 , the integrals F 1 
and G-y vanish in the far-field and amount to very small values when evaluated 
numerically around the core. Thus, we conclude that to leading order in the 
non-potential terms (and the wavenumber q ) the velocity of an individual 
PHD is well approximated by 



L xx ± xy 



\ lyx lyy J 



G n 



+ 



G, 



a 1 F 1 + a 2 F 2 
a 3 G 3 



(24) 



To evaluate the mobility tensor on the left-hand side of (24) it is not sufficient 
to insert the solution for the stationary defect in (19) since it leads to integrals 
that diverge in the far field. To regularize this singularity the solution for the 
moving defect has to be used, at least in the far field [34,32,31]. To leading 
order, the non-potential terms can be neglected in the mobility tensor. Then 
its off-diagonal terms are much smaller than the diagonal terms [31]. Using 
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Fig. 14. Defect velocity due to ot\ t 2 (glide) and due to (climb), respectively, for 
q = 0, ft = 1, 7 = 0.2, (3 = 0, and v = 2. 

the numerically determined defect solution of the full equations (1,2), we also 
find I xy to be much smaller than either I xx or I yy (by a factor of 100) when we 
only compute the integral within a box enclosing the defect core and neglect 
the far-field contribution. 

The effect of the non-linear gradient terms can now be summarized as follows: 
since the off-diagonal terms in the mobility tensor are small, the contributions 
of cti and a.2 to the velocity of a PHD with dislocations in amplitudes A2 and 
A3 are in the rii-direction ("glide"), while a% causes such a PHD to travel 
in the fi-direction ("climb"). Furthermore, like the mean flow the nonlinear 
gradient terms contribute to a shift of the wavenumber at which the defect is 
stationary away from the band center q = 0. 

The above conclusion is confirmed in direct numerical simulations of (1,2). In 
Fig. 14 we plot the velocity of a PHD with charges (0, — , +) for fi = 1, q = 0, 
7 = 0.2, and (3 = as a function of the strength of the nonlinear gradient 
terms. As expected from the discussion of (24), if only oti or «2 are non-zero 
the defect glides, whereas it climbs if only a 3 is present. The velocities change 
sign if the charges of the PHD are reversed. 

Fig. 14 indicates that the linear dependence of the velocity on the is re- 
stricted to a range \oti\ < 0.2. It is found, however, that the direction of the 
defect motion remains the same for aj outside that range, i.e. a.\ and 0:2 cause 
defects to glide while defects climb due to a^. In comparison, in previous nu- 
merical simulations we found that the mean flow causes the PHD's to perform 
a combined climb-glide motion [14]. This indicates that both F@ and Gp are 
non-zero in equation (24). Fig. 15 shows the dependence of the defect velocity 
on the wavenumber q. For a 3 7^ 0, the defect performs a mixed climb-glide 
motion if the wave number is not at the band center, i.e. if q 7^ 0. 
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Fig. 15. Dependence of the defect velocity on the wavenumber q for \i = 1, v = 2 
and f3 = 0. cxi = except as indicated: a\ = 0.3 (circles), = 0.3 (triangles), 
03 = 0.3 (squares). Dotted line is for the variational case a\ = a<i = 03 = 0. 



6 Conclusion 



The variational character of the minimal Ginzburg-Landau equations for steady 
hexagon patterns can be broken in various ways. Quite generally, at cubic or- 
der two nonlinear gradient terms arise [35,36]. If the chiral symmetry of the 
system is broken (e.g. by rotation) a third nonlinear gradient term is possible 
[8]. In Rayleigh-Benard and in Marangoni convection the mean flow driven by 
deformations of the pattern introduce a further non- variational term in the 
form of a non-local coupling of the roll modes (e.g. [16,37]). In this paper we 
have studied the combined effect of the mean flow and of rotation (includ- 
ing the respective nonlinear gradient term) on the stability and dynamics of 
hexagon patterns as well as on the stability and motion of their penta-hepta 
defects. 

Rotation induces a coupling of the two long-wave phase modes that can gen- 
erate a long-wave oscillatory instability, which in the absence of nonlinear 
gradient terms is, however, usually screened by a steady short-wave instabil- 
ity [8]. The mean flow is only driven by the transverse phase mode [14]. Our 
results indicate that it suppresses the rot at ion- induced coupling of the phase 
modes that leads to oscillatory behavior. 

The various non- variational terms influence the motion of penta-hepta defects 
in different ways. While the nonlinear gradient terms that preserve the chiral 
symmetry contribute to a gliding motion of the defects, the nonlinear gra- 
dient term introduced by rotation induces a climbing motion. This is to be 
contrasted with the effect of the mean flow, which leads to a mixed climbing- 
gliding motion. 

The most interesting result of this paper is associated with the fact that the 
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nonlinear gradient terms lead to a separation of the two dislocations mak- 
ing up a penta-hepta defect and can destabilize it [22]. More specifically, in 
simulations we find that penta-hepta defects induce the nucleation of new dis- 
locations in the defect-free amplitude. Even a weak mean flow can enhance this 
instability significantly. Moreover, with a sufficiently strong nonlinear gradient 
term that breaks the chiral symmetry (0:3) the induced nucleation can lead 
to an ever-increasing transverse wavevector component of the hexagon pat- 
terns. As in the case without mean flow [8], we interpret this growth, which 
eventually leads to the break-down of the Ginzburg-Landau equations, as the 
signature of a persistent precession of the disordered pattern on average. In 
Fourier space, the wave-vector spectrum of such a precessing pattern would 
drift along the critical circle. In the lowest-order Ginzburg-Landau equations 
used here the critical circle is replaced by its tangents at each of the three 
modes making up the hexagons, which are transverse to the respective wave- 
vectors. We expect therefore that in this regime the physical system would 
exhibit persistent penta-hepta defect chaos driven by induced nucleation. 

To overcome the limitations of the Ginzburg-Landau equations, we are cur- 
rently investigating penta-hepta defect chaos using a suitably modified Swift- 
Hohenberg-type equations coupled to a mean flow [27]. As expected from the 
simulations of the Ginzburg-Landau equations presented in the present paper, 
the penta-hepta defect chaos persists and on average the disordered hexagon 
patterns precess indefinitely. As in the Ginzburg-Landau case, the chaotic 
state is due to the induced nucleation of penta-hepta defects, which is ap- 
parently only possible when the nonlinear gradient terms are included. These 
simulations also indicate that the induced nucleation is not contingent on the 
oscillatory sideband instability; rather it is the separation of the dislocations in 
a penta-hepta defect that renders them unstable and induces the nucleation. 
To obtain persistent chaos the resulting defect proliferation has to be strong 
enough compared to the annihilation rate and apparently the chiral symmetry 
has to be broken. 
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